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Abstract 

We study the deep interplay between geometry of quadrics in d-dimen- 
sional space and the dynamics of related integrable billiard systems. Var¬ 
ious generalizations of Poncelet theorem are reviewed. The corresponding 
analytic conditions of Cayley’s type are derived giving the full description 
of periodical billiard trajectories; among other cases, we consider billiards 
in arbitrary dimension d with the boundary consisting of arbitrary num¬ 
ber k of confocal quadrics. Several important examples are presented in 
full details demonstrating the effectiveness of the obtained results. We 
give a thorough analysis of classical ideas and results of Darboux and 
methodology of Lebesgue, and prove their natural generalizations, ob¬ 
taining new interesting properties of pencils of quadrics. At the same 
time, we show essential connections between these classical ideas and the 
modern algebro-geometric approach in the integrable systems theory. 
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1 Introduction 


In his Traits des proprietes projectives des figures [34], Poncelet proved one 
of most beautiful and most important claims of the 19th century geometry. 
Suppose that two ellipses are given in the plane, together with a closed polygonal 
line inscribed in one of them and circumscribed about the other one. Then, 
Poncelet theorem states that infinitely many such closed polygonal lines exist 
- every point of the first ellipse is a vertex of such a polygon. Besides, all 
these polygons have the same number of sides. Poncelet’s proof was purely 
geometrical, synthetic. Later, using the addition theorem for elliptic functions, 
Jacobi gave another proof of the theorem m- Essentially, Poncelet theorem is 
equivalent to the addition theorem and Poncelet’s proof represents a synthetic 
way of deriving the group structure on an elliptic curve. Another proof, in 
a modern, algebro-geometrical manner, can be found in Griffiths’ and Harris’ 
paper [^ . There, they also gave an interesting generalization of the Poncelet 
theorem to the three-dimensional case, considering polyhedral surfaces both 
inscribed and circumscribed about two quadrics. 

A natural question connected with Poncelet theorem is to find an analytical 
condition determining, for two given conics, if an n-polygon inscribed in one 
and circumscribed about the second conic exists. In a short paper Uni, Cay¬ 
ley derived such a condition, using the theory of Abelian integrals. Inspired 
by this paper, Lebesgue translated Cayley’s proof to the language of geome¬ 
try. Lebesgue’s proof of Cayley’s condition, derived by methods of projective 
geometry and algebra, can be found in his book Les coniques ED. Griffiths and 
Harris derived Cayley theorem by finding an analytical condition for points of 
finite order on an elliptic curve |24] . 

It is worth emphasizing that Poncelet, in fact, proved a statement that is 
much more general than the famous Poncelet theorem EllMj, then deriving 
the latter as a corollary. Namely, he considered n -|- 1 conics of a pencil in the 
projective plane. If there exists an n-polygon with vertices lying on the first of 
these conics and each side touching one of the other n conics, then infinitely 
many such polygons exist. We shall refer to this statement as Complete Poncelet 
theorem (CPT) and call such polygons Poncelet polygons. We are going to follow 
here mostly the presentation of Lebesgue from m, which is, as we learned from 
M. Berger, quite close to one of two Poncelet’s original proofs. 

A nice historical overview of the Poncelet theorem, together with modern 
proofs and remarks is given in |9| . Various classical theorems of Poncelet type 
with short modern proofs reviewed in , while the algebro-geometrical approach 
to families of Poncelet polygons via modular curves is given in ISIEH]. 

Poncelet theorem has a nice mechanical interpretation. Elliptical billiard 
|30j is a dynamical system where a material point of the unit mass is moving 
with a constant velocity inside an ellipse and obeying the reflection law at the 
boundary, i.e. having congruent impact and reflection angles with the tangent 
line to the ellipse at any bouncing point. It is also assumed that the reflection is 
absolutely elastic. It is well known that any segment of a given elliptical billiard 
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trajectory is tangent to the same conic, confocal with the boundary m- If 
a trajectory becomes closed after n reflections, then Poncelet theorem implies 
that any trajectory of the billiard system, which shares the same caustic curve, 
is also periodic with the period n. 

Complete Poncelet theorem also has a mechanical meaning. The configura¬ 
tion dual to a pencil of conics in the plane is a family of confocal second order 
curves [4]. Let us consider the following, a little bit unusual billiard. Suppose 
n confocal conics are given. A particle is bouncing on each of these n conics 
respectively. Any segment of such a trajectory is tangent to the same conic 
confocal with the given n curves. If the motion becomes closed after n reflec¬ 
tions, then, by Complete Poncelet theorem, any such a trajectory with the same 
caustic is also closed. 

The statement dual to Complete Poncelet theorem can be generalized to 
the d-dimensional space Suppose vertices of the polygon xiX 2 ■ ■ ■ Xn are 
respectively placed on confocal quadric hyper-surfaces Qi, Q2, • ■ •, Qn in the d- 
dimensional Eucledean space, with consecutive sides obeying the reflection law 
at the corresponding hyper-surface. Then all sides are tangent to some quadrics 
Q^, • ■ ■, Q'^~^ confocal with {Qi}; for the hyper-surfaces {Qi,Q^}, an infinite 
family of polygons with the same properties exist. 

But, more than one century before these quite recent results, Darboux proved 
the generalization of Poncelet theorem for a billiard within an ellipsoid in the 
three-dimensional space m- It seems that his work on this topic is completely 
forgot nowadays. 

It is natural to search for a Cayley-type condition related to some of gen¬ 
eralizations of Poncelet theorem. The authors derived such conditions for the 
billiard system inside an ellipsoid in the Eucledean space of arbitrary finite 
dimension [H HO]. In our recent note im, algebro-geometric conditions for 
existence of periodical billiard trajectories within k quadrics in d-dimensional 
Euclidean space were announced. The aim of the present paper is to give full ex¬ 
planations and proofs of these results together with several important examples 
and improvements. The second important goal of this paper is to offer a thor¬ 
ough historical overview of the subject with a special attention on the detailed 
analysis of ideas and contributions of Darboux and Lebesgue. While Lebesgue’s 
work on this subject has been, although rarely, mentioned by experts, on the 
other hand, it seems to us that relevant Darboux’s ideas are practically un¬ 
known in the contemporary mathematics. We give natural higher dimensional 
generalizations of the ideas and results of Darboux and Lebesgue, providing the 
proofs also in the low-dimensional cases if they were omitted in the original 
works. Beside other results, interesting new properties of pencils of quadrics are 
established - see Theorems [5] and [TUI The latter gives a nontrivial generalization 
of the Basic Lemma. 

This paper is organized as follows. In the next section, a short review of 
Lebesgue’s results from m is given, followed by their application to the case 
of the billiard system between two confocal ellipses. Section [3] contains algebro- 
geometric discussions which will be applied in the rest of the paper. In Section 
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01 we give analytic conditions for periodicity of billiard motion inside a do¬ 
main bounded by several confocal quadrics in the Euclidean space of arbitrary 
dimension. The complexity of the problem of billiard motion within several 
quadrics is well known, even in the real case, and it is induced by multivalued¬ 
ness of the billiard mapping. Thus, to establish a correct setting of the problem, 
we introduce basic notions of reflections from inside and from outside a quadric 
hyper-surface, and we define the billiard ordered game. The corresponding close¬ 
ness conditions are derived, together with examples and discussions. In Section 
01 we consider the elliptical billiard as a discrete-time dynamical system, and, 
applying the Veselov’s and Moser’s algebro-geometric integration procedure, we 
derive the periodicity conditions. The obtained results are compared with those 
from Section 01 In Section [6l we give an algebro-geometric description of pe¬ 
riodical trajectories of the billiard motion on quadric hyper-surfaces, we study 
the behaviour of geodesic lines after the reflection at a confocal quadric and 
derive a new porism of Poncelet type. In Section 0 we define the virtual reflec¬ 
tion configuration, prove Darboux’s statement on virtual billiard trajectories, 
generalize it to arbitrary dimension and study related geometric questions. In 
Section [8l we formulate and prove highly nontrivial generalization of the Basic 
Lemma (Lemma 01), giving a new important geometric property of dual pencils 
of quadrics. In that section, we also introduce and study the generalized Cayley 
curve, a natural higher-dimensional generalization of the Cayley cubic studied 
by Lebesgue. In this way, in Section 01 the most important tools of Lebesgue’s 
study are generalized. Further development of this line will be presented in 
separate publication [22j . In Appendix I, we review some known classes of in- 
tegrable potential perturbations of elliptical billiards, emphasizing connections 
with Appell hypergeometric functions and Liouville surfaces. Finally, in Ap¬ 
pendix 2, we present the related Darboux’s results considering a generalization 
of Poncelet theorem to Liouville surfaces, giving a good basis for a study of 
the geometry of periodic trajectories appearing in the perturbed systems from 
Appendix 1. 


2 Planar Case: d = 2, k — Arbitrary 

First of all, we consider the billiard system within k confocal ellipses in the 2- 
dimensional plane. In such a system, the billiard particle bounces sequentially 
of these confocal ellipses. We wish to get the analytical description of periodical 
trajectories of such a system. 

Following Lebesgue, let us consider polygons inscribed in a conic F, whose 
sides are tangent to ri,...,Ffc, where F,ri,...,Ffe all belong to a pencil of 
conics. In the dual plane, such polygons correspond to billiard trajectories 
having caustic F* with bounces on FJ,..., F^. The main object of Lebesgue’s 
analysis is the cubic Cayley curve, which parametrizes contact points of tangents 
drawn from a given point to all conics of the pencil. 
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2.1 Full Poncelet Theorem 


Basic Lemma. Next lemma is the main step in the proof of full Poncelet 
theorem. If one Poncelet polygon is given, this lemma enables us to construct 
every Poncelet polygon with given initial conditions. Also, the lemma is used 
in deriving of a geometric condition for the existence of a Poncelet polygon. 

Lemma 1 [31] Let T he a pencil of conics in the projective plane and T a conic 
from this pencil. Then there exist quadrangles whose vertices A, B, C, D are on 
P such that three pairs of its non-adjacent sides AB, CD; AC, BD; AD, BC 
are tangent to three conics of T. Moreover, the six contact points all lie on a 
line A. Any such a quadrangle is determined by two sides and the corresponding 
contact points. 

Let P, Pi, r 2 , Pa be conics of a pencil and ABC a Poncelet triangle corre¬ 
sponding to these conics, such that its vertices lie on P and sides AB, BC, CA 
touch Pi, P 2 , P 3 respectively. This lemma gives us a possibility to construct 
triangle ABD inscribed in P whose sides AB, BD, DA touch conics Pi, P 3 , P 2 
respectively. In a similar fashion, for a given Poncelet polygon, we can, applying 
Lemma 1, construct another polygon which corresponds to the same conics, but 
its sides are tangent to them in different order. 

Circumscribed and Tangent Polygons. Let a triangle ABC be inscribed 
in a conic P and sides BC, AC, AB touch conics Pi, P 2 , P 3 of the pencil B at 
points M, N, P respectively. According to the Lemma, there are two possible 
cases: either points M, N, P are collinear, when we will say that the triangle 
is tangent to Pi, P 2 , P 3 ; or the line MN is intersecting AB at a point S which 
is a harmonic conjugate to P with respect to the pair A, B, then we say that 
the triangle is circumscribed about Pi, P 2 , P 3 . 

Let ABCD ... KL he a polygon inscribed in P whose sides touch conics Pi, 

..., P„ of the pencil F respectively. Denote by {AC) a conic such that AABC 
is circumscribed about Pi, P 2 , {AC), by {AD) a conic such that AACD is 
circumscribed about {AC), P 3 , {AD). Similarly, we find conics {AE), ..., 
{AK). The triangle AKL can be tangent to the conics {AK), Pn-i, Pn or 
circumscribed about them, and we will say that ABCD ... KL is tangent or, 
respectively, circumscribed about conics Pi, ..., P„. 

Further, we will be interested only in circumscribed polygons. The Poncelet 
theorem does not hold for tangent triangles nor, hence, for tangent polygons 
with greater number of vertices. 

Theorem 1 (Complete Poncelet theorem) Let conics P, Pi, ..., P„ belong to 
a pencil F. Lf a polygon inscribed in P and circumscribed about Pi, ..., P„ 
exists, then infinitely many such polygons exist. 

To determine such a polygon, it is possible to give arbitrarily: 

1) the order which its sides touch Pi, ..., P„ in; let the order be: P'l, ..., Pj^; 

2 ) a tangent to Pj containing one side of the polygon; 
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3) the intersecting point of this tangent with T which will belong to the side 
tangent to 

The proof is given in m- 

2.2 Cayley’s Condition 

Representation of Conics of a Pencil by Points on a Cubic Curve. Let 

pencil T of conics be determined by the curves (7 = 0 and F = 0. The equation 
of an arbitrary conic of the pencil is (7 + AF = 0. 

Let P + An = 0 be the equation of the corresponding polar lines from the 
point n e F. The geometric place of contact points of tangents from A with 
conics of the pencil is the cubic C : CA — TP = 0. On this cubic, any conic 
of P is represented by two contact points, which we will call representative 
points of the conic. The line determined by these two points passes through 
point Z : P = 0, H = 0. There exist exactly four conics of the pencil whose 
representative points coincide: the conic F and three degenerate conics with 
representative points d, a, f3, 7 . Lines ZA, Za^ Zfd, Z 7 are tangents to C 
constructed from Z. The tangent line to cubic C at point Z is a polar of point 
A with respect to the conic of the pencil which contains Z. 

Condition for Existence of a Poncelet Triangle. If triangles inscribed in 
F and circumscribed about Fi, F 2 , F 3 exist, we will say that the conics Fi, F 2 , 
F 3 are joined to F. In this case, CPT states there is six such triangles with the 
vertex d e F. Let ABC be one of them. Side AB, denote it by 1, touches Fi 
in point mi. Also, it touches another conic of the pencil, denote it by (/), in 
point M. Side AC, denote it by 2’, touches F 2 in m^. Consider the quadrangle 
ABCD determined by AB, AC and contact points M, m^. Line Mm 2 meets 
BC at p 3 , its point of tangency to F3, and meets AD (which we will denote by 
S’) at the point m3 of tangency to F3. Triangle ABD is circumscribed about 
Pi, F 2 , Fa- 

Similar ly, triangle ACE can be obtained by construction of quadrangle 
ABCE determined by AB, BC and contact points mi, Line AE touches 
F 3 at point m 3 S mi^ 3 . Denote this line by 3. Triangles with sides S’, 2 and 3, 
1 ’ are constructed analogously. 

There is exactly six tangents from A to conics Fi, F2, F3. We have divided 
these six lines into two groups: 1,^,5 i l’,2’,3’. Two tangents enumerated by 
different numbers and do not belong to the same group, determine a Poncelet 
triangle. 

Cubic C and the cubic consisting of lines m^M, m2m2, m^m'^ have simple 
common points mi, M, m 2 , m' 2 , m3, m3. A, and point Z as a double one. 
A pencil determined by these two cubics contains a curve that passes through 
a given point of line Mm 2 , different from M, m 2 , m3. This cubic has four 
common points with line Mm' 2 , so it decomposes into the line and a conic. 
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Thus, mi, m 2 , m 3 are intersection points, different from A and Z, of a conic 
which contains A and touches cubic C at Z. 

Converse also holds. 

Let an arbitrary conic that contains point A and touches cubic C at Z be 
given. Denote by mi, m 2 , m 3 remaining intersection points of the curve C with 
this conic. Each of the lines miZ, m 2 Z, m^Z has another common point with 
the cubic C; denote them by m'l, m^, m 3 respectively. By definition of the curve 
C, we have that mi,m[; m 2 ,m 2 ; m 3 ,m 3 are pairs of representative points of 
some conics Li, r 2 , r 3 from the pencil T. Line Am-i, besides being tangent to 
Li at mi, has to touch another conic from the pencil T. Take that it is tangent 
to a conic (/) at M. 

Now, in a similar fashion as before, we can conclude that points M, m^, m 3 
are collinear. Applying Lemma 1, it is easily deduced that conics Li, r 2 , r 3 are 
joined to L. 

So, we have shown the following: systems of three joined conics are deter¬ 
mined by systems of three intersecting points of cubic C with conics that contain 
point A and touch the curve C at Z. 

Cayley’s Cubic. Let D{X) be the discriminant of conic C + AL = 0. We will 
call the curve 

Co : T" = D{X) 

Cayley’s cubic. Representative points of conic C + AT = 0 on Cayley’s cubic 
are two points that correspond to the value X = A. 

The polar conic of the point Z with respect to cubic C passes through the 
contact points of the tangents ZA, Za, Z/3, Z 7 from Z to C. Thus, points a, /3 ,7 
are representative points of three joined conics from the pencil IF. Those three 
conics are obviously the decomposable ones. Corresponding values A diminish 
D{X), and these three representative points on Cayley’s cubic Cq lie on the line 

y = 0. 

Using the Sylvester’s theory of residues, we will show the following: 

Let three representative points of three conics of pencil T be given on the 
Cayley’s cubic Cq. Condition for these conics to be joined to the conic T is that 
their representative points are collinear. 

Sylvester’s Theory of Residues. When considering algebraic curves of genus 
1 , like the Cayley’s cubic is here, Abel’s theorem can always be replaced by 
application of this theory. 

Proposition 1 Let a given cubic and an algebraic curve of degree m-\- n meet 
at 3(m + n) points. If there is 3m points among them which are placed on a 
curve of degree m, then the remaining 3n points are placed on a curve of degree 


n. 



If the union of two systems of points is the complete intersection of a given 
cubic and some algebraic curve, then we will say that these two systems are 
residual to each other. Now, the following holds: 

Proposition 2 If systems A and A' of points on a given cubic curve have a 
common residual system, then they share all residual systems. 

Proof. Suppose ,8 is a system residual to both A, A' and B' is residual to A. 
Then AU A' is residual to BUB' , i.e. the system A^J A' U A^J A' is a complete 
intersection of the cubic with an algebraic curve. Since ^ U 8 is also such an 
intersection, it follows, by the previous proposition, that A' and B' are residual 
to each other. □ 

Let us note that this proposition can be derived as a consequence of Abel’s 
theorem, for a plane algebraic curve of arbitrary degree. However, if the degree 
is equal to three, i.e. the curve is elliptic, Proposition 2 is equivalent to Abel’s 
theorem. 

Condition for Existence of a Poncelet Polygon. Let conics L, Li,..., r„ 

be from a pencil. If there exists a polygon inscribed in P and circumscribed 
about Pi,... , r„, we are going to say that conics Ti,..., r„ are joined to P. 
Then, similarly as in the case of the triangle, it can be proved that tangents 
from the point A € P to Pi,..., P„ can be divided into two groups such that any 
Poncelet n-polygon with vertex A has exactly one side in each of the groups. 

This division of tangents gives a division of characteristic points of conics 
Pi, ..., P„ into two groups on C and, therefore, a division into two groups on 
Cayley’s cubic Cq: .. .and _ 

Let ABCD .. .KL he a, Poncelet polygon, and let (AC), {AD), ... be conics 
determined like in the definition of a circumscribed polygon. Let c, d,5', ... be 
a corresponding characteristic points on Cq, such that triples 1, 2, c; 7 , 3, d; 6, 
4, e are characteristic points of the same group with respect to corresponding 
conics. 

Points 1, 2, c are collinear, as 7 , 3, d are. Thus, 1, 2, 3, d are residual 
with c, 7 . Line 07 contains point Z, so system c, 7 is residual with Z, too. It 
is possible to show that Z is a triple point of curve Cq and it follows that it is 
residual with system Z, Z. This implies that points 1,2,3, d, Z, Z are placed on 
a conic. 

If we take a coordinate system such that the tangent line to Co at Z is the 
infinite line and the axis Oy is line AZ, we will have: 

four conics are joined to P if and only if their characteristic points of the 
same group are on a parabola with the asymptotic direction Oy. 

Continuing deduction in the same manner, we can conclude: 3n — p points of 
the cubic Co are characteristic points of same group for in—p conics (1 < p < 3) 
joined to P, if and only if these points are placed on a curve of degree n which 
has Oy as an asymptotic line of the order p. 
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Cayley’s Condition. Let Cq : = D{x) be the Cayley’s cubic, where D{x) is 

the discriminant of the conic C + a;r = 0 from pencil T. A system of n conics 
joined to F is determined by n values x if and only if these n values are abscissae 
of intersecting points of Co and some algebraic curve. Plugging D(x) instead of 
in the equation of this curve, we obtain: 

P{x)y + Q{x) = 0, 


that is _ 

P{x)\JD{x) + Q{x) = 0. 

From there: 

D{x){aoxP~‘^ + aixP~^ H-h 0 ^- 2 ) + (box^ + bix^~^ H- \-bp) =0, n = 2p; 

a/ D{x){aox^~^ +aix^~‘^ + ■ ■ ■+ap-i) + {box^+bix^~^ + - ■ ■+bp) =0, n = 2p+l. 

If Ai,..., Afc denote parameters corresponding to Fi,..., Ffc respectively, 
then existence of a Poncelet polygon inscribed in F and circumscribed about 
Fi,..., Ffc is equivalent to: 


1 Ai XI 

1 Afe XI 

for k = 2p; 

1 Ai Af 

1 Afe XI 

for fc = 2 p + 1 . 


A? 




A? 


K 


\/^(Ai) 

Xi^D(X^) .. 

. Arv^(Ai) 

^jD{Xk) 

Xk\jD{Xk) .. 

. Xr^y/D{Xk) 

\/^(Ai) 

XisjD{X^) .. 

. Arv^(Ai) 

^jD{Xk) 

Xk\jD{Xk) .. 

. Xr^y/D{Xk) 


There exists an n-polygon inscribed in F and circumscribed about C if and 
only if it is possible to find coefficients oq, ui, •..; bo, 61 , ... such that function 
P{x)\JD{x) + Q{x) has X = 0 as a root of the multiplicity n. 

For n = 2p, this is equivalent to the existence of a non-trivial solution of the 
following system: 


aoCs 

+ 

aiC4 

+ .. 

. + 


= 0 

aoCi 

+ 

ctiCo 

+ .. 

. + 

(lp—2Cp-\.2 

= 0 

aoCp+i 

+ 

aiC'p+2 

+ .. 

. + 

ap-2C2p-l 

= 0 , 


where 

a/ D{x) = A + Bx + C2X^ + Cox^ + - 
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Finally, for n = 2p, we obtain the Cayley’s condition 


Similarly, for n 


C 3 

C 4 

■ Cp+i 


C 4 

C 5 

■ Cp+2 

= 0 , 

Cp+i 

Cp+2 ■ • 

C 2 P -1 


1 , we obtain: 



C 2 

C 3 

■ ■ Cp+i 


M 

C 4 

■ ■ Cp+2 

= 0 . 

Cp+i 

Cp+2 ■ 

.. C2p 



These results can be directly applied to the billiard system within an el¬ 
lipse: to determine whether a billiard trajectory with a given confocal caustic 
is periodic, we need to consider the pencil determined by the boundary and the 
caustic curve. 


2.3 Some Applications of Lebesgue’s Results 

Now, we are going to apply the Lebesgue’s results to billiards systems within 
several confocal conics in the plane. 

Consider the dual plane. The case with two ellipses, when the billiard tra¬ 
jectory is placed between them and particle bounces to one and another of them 
alternately, is of special interest. 


Corollary 1 The condition for the existence of 2m-periodic billiard trajectory 
which bounces exactly m times to the ellipse Tj; = C* and m times to F^ = 
(C -I- yF)*, having F* for the caustic, is: 


/ /o(0) /i(0) 

m) fm 


/2m-l(0) \ 

/2m-l(0) 1 


det 


foil) Ml) 

foil) Ml) 


ftl'm 

/2m-l(7) 

fM-iii) 


\ft~'\i) f^^-^\i) 


Am-1) 

J2m-1 


ii) 


J 


where fj = , (0 < j < m), fm+i = ^ \/ DM), (1 < i < m — 1). 


We consider a simple example with four bounces on each of the two conics. 


Example 1 The condition on a billiard trajectory placed between ellipses F* 
and r 2 , to be closed after 4 alternate bounces to each of them is: 

det X = 0 , 
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where the elements of the 3x3 matrix X are: 

Xii = -4Bo + Bij + 4Co + 3(717 + 2(^27^ + <737^ 

Xi 2 = -3Bo + B^-f + 3(7o + 2 Ci7 + C 27 " 

Xi 3 = -2So + Si 7 + 2(7o + Ci7 

X 21 = - 6 S 0 + 527 ^ + 6(7o + 6 ( 7 i 7 + 4(727^ + 3(737^ 

X 22 = - 6 S 0 + Si 7 + 527^ + 6(7o + 4(7 i7 + 3(727^ 

X 23 = -5Bo + 2 Bi 7 + ^ 27 ^ + 5(7o + 3(7i7 

X31 = -4Bo + S373 + 4(70 + 4(7 i7 + 4(727^ + 3(737^ 

X 32 = -4Bo + ^ 27 ^ + B 373 + 4(7o + 4(7 i7 + 3(727^ 

X 33 = —4_Bo + Bi^ + - 627 ^ + + 4(7o + 3 ( 7 i 7 , 

with (7i, Bi being coefficients in the Taylor expansions around x = 0 and a; = 7 
respectively: 

A^/ D[x) = Cq + Cix + C 2 x'^ + ..., 

\/B>{x) = Bq + Bi{x - j) + B 2 {x - + - 

On Figure [H we see a Poncelet octagon inscribed in F and circumscribed 
about Fi and F 2 . In the dual plane, the billiard trajectory that corresponds to 
this octagon, has the dual conic F* as the caustic (see Figure [2]). 




Figure 1: A Poncelet octagon 
whose sides touch the conics Fi 
and F 2 alternately 


Figure 2: A closed billiard trajec¬ 
tory in the domain bounded by 
two confocal ellipses 


3 Points of Finite Order on the Jacobian of a 
Hyperelliptic Curve 

In order to prepare the algebro-geometric background for the rest of the article, 
in this section we are going to give the analytical characterization of some classes 
of finite order divisors on a hyperelliptic curve. 
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Let the curve C be given by 

= {x - xi).. .{x - X 2 g+i), Xi^Xj when i ^ j. 

It is a regular hyperelliptic curve of genus g, embedded in Let J{C) be its 
Jacobian variety and 

A : C^J{C) 

the Abel-Jacobi map. 

Take E to be the point which corresponds to the value x = oo, and choose 
A{E) to be the neutral in J{C). According to the Abel’s theorem A{Pi) + 
■ ■ ■ + A{Pn) = 0 if and only if there exists a meromorphic function / with zeroes 
Pi,... ,Pn and a pole of order n at the point E. Let C{nE) be the vector space 
of meromorphic functions on C with a unique pole E of order at most n, and 
fi,..., fk a basis of C(nE). The mapping 

F-.C^P>^-\ X^[h{X),...,h{X)] 

is a projective embedding whose image is a smooth algebraic curve of degree 
n. Hyperplane sections of this curve are zeroes of functions from C{nE). Thus, 
the equality nA{P) = 0 is equivalent to: 


/ fl{p) 

f2{P) 

■ fk{p) \ 



fiip) 

mp) 

■ fkip) 

< k. 

(1) 

V (P) 

ft~'\p) ■ 

. ft~'\p)) 




Lemma 2 For n < 2g, there does not exist a point P on the curve C, such that 
nA{P) = 0 and P ^ E. 

Proof. Let P be a point on C, P ^ E, and x = xq its corresponding value. 
Consider the case of n even. Since P is a branch point of a hyperelliptic curve, its 
Weierstrass gap sequence is 1, 3, 5,..., 2g — 1 Now, applying the Riemann- 
Roch theorem, we obtain dim£(nP) = n/2 + 1. Choosing a basis l,x,..., 
for C(nE), and substituting in ([T]), we come to a contradiction. □ 


Lemma 3 Let P{xo,yo) he a non-branching point on the curve C. For n > 2g, 
equality nA{P) = 0 is equivalent to: 




^ -^m+1 


Bm 

■ ■ ■ Bg+2 

\ 

rank 

Bm+2 


Bm-\-l 

■ ■ ■ Bg + 3 




\ B2m- 

1 

B2m-2 

■ ■ ■ BjYi-\-g j 


/ 

^m+1 

Bm • • • 

Bg + l \ 


rank 


Bm+2 

Bm-\-l • ■ ■ 

Bg+2 

< 


1 

^2m 


Bm-\-g j 


and yj [x 

- Xi). .. 

{x 

- X2g+l) 

= Bo-\-Bi( 

X — 


< m — g, when n = 2m, (2) 


< m — g 1, when n = 2m 1, 
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Proof. This claim follows from previous results by choosing a basis for C{nE): 

l,x,..., a;™, y, xy,..., x^~^~^y if n = 2 m, 

1 , X,, x™, y,xy,.x'^~^y if n = 2 m + 1 , 

similarly as in [^. □ 

In the next lemma, we are going to consider the case when the curve C is 
singular, i.e. when some of the values Xi, X 2 , ■ ■ ■, 2 ^ 29+1 coincide. 

Lemma 4 Let the curve C be given by 

y'^ = {X - Xi) . . .{X - X2g+l), Xi- X 2 . 2 : 29+1 ^ 0, 

Pq one of the points corresponding to the value x = Q and E the infinite point 
on C. Then 2nPo ~ 2nE is equivalent to m, where 

y = {x - xi).. .{x - X 2 g+i) = BqP Bix + B 2 X^ + ... 
is the Taylor expansion around the point Pq . 

Proof. Suppose that, among 2 : 1 ,..., 2 : 29 + 1 , only X 2 g and 2 : 29+1 have same values. 
Then ( 2 : 29 , 0) is an ordinary double point on C. The normalization of the curve 
C is the pair (C, tt), where C is the curve given by: 

C : y"^ = {x - xi)... {x - 2:29-1), 

and TT : C C is the projection: 

(i, y)iBU {x = x, y={x- X 2 g)y). 

The genus of C is (jf — 1. 

Denote by A and B points on C which are mapped to the singular point 
(x 2 g, 0) G C by the projection tt. Any other point on C is the image of a unique 
point of the curve C. Let 


tt{E) = E, tt{Po) = Po. 

The relation 2nPo ^ 2nE holds if and only if there exists a meromorphic 
function / on C, f € C{2nE), having a zero of order 2n at Pq and satisfying 

fiA) = fiB). 

For n < g — 1, according to Lemma [31 2nE ^ 2nPo cannot hold. For n> g, 
choose the following basis of the space C{2nE): 

® ■ ■ ■ j fn—g—1 ^ 

where 1, /i,..., fn-g-i is a basis of C{2nE) as in the proof of LemmajSl 
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Since y is the only function in the basis which has different values at points 
A and we obtain that the condition 

2nE ~ 2npQ 

is equivalent to dH). 

Cases when C has more singular points or singularities of higher order, can 
be discussed in the same way. □ 


Lemma 5 Let the curve C he given by 

y"^ = {X - Xi) ... {X - X2g+2), 


with all Xi distinct from 0, and Q+, Q- the two points on C over the point 
X = 0. Then nQ+ ^ nQ- is equivalent to: 


^ Bg + 2 -Bff+S 
Bg+S Bg+4 


Bn+l ^ 
Bn+2 


rank 


< n — g and n > g, (3) 


\ Bg4.n 


B2n-1 j 


where y = y/{x — xi)... {x — X 2 g+ 2 ) = Bq + BiX + B 2 X^ + ... is the Taylor 
expansion around the point Q-. 


Proof. C is a hyperelliptic curve of genus g. The relation n(5+ = nQ- means 
that there exists a meromorphic function on C with a pole of order n at the 
point Q+, a zero of the same order at Q- and neither other zeros nor poles. 
Denote by LinQ^.) the vector space of meromorphic functions on C with a 
unique pole Q+ of order at most n. Since Q+ is not a branching point on the 
curve, dimL(n(5+) = 1 for n < g, and dim L{nQ+) = n — g + 1, for n > g. 
In the case n < g, the space contains only constant functions, and the 

divisors nQ+ and nQ- can not be equivalent. If n > 5 + 1, we choose the 
following basis for L{nQ+): 

1) /li ■ • • I fn—g, 


where 

r _y- Bq- BiX - 

~ a;S+fe 

Thus, n(5+ = nQ- if there is a function / S L{nQ+) with a zero of order n 
at (5-, i.e., if there exist constants ao,..., an-g, not all equal to 0, such that: 


ao + aifi{Q-) + 

o^ifiiQ-) + 

^\Q-) + 


T ^n—gfn—g{Q—') — 0 

+ an-gfn-glQ-) = 0 

+ ar.-gft~g^\Q-) = 0 . 
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Existence of a non-trivial solution to this system of linear equations is equivalent 
to the condition ([3]). 

When some of the values xi,... ,X 2 g +2 coincide, the curve C is singular. 
This case can be considered by the procedure of normalization of the curve, as 
in Lemma 4. The condition for the equivalence of the divisors nQ^ i nQ-, in 
the case when C is singular, is again □ 

4 Periodic Billiard Trajectories inside k Confo- 
cal Quadrics in 


Darboux was the first who considered a higher-dimensional generalization of 
Poncelet theorem. Namely, he investigated light-rays in the three-dimensional 
case {d = 3) and announced the corresponding complete Poncelet theorem in 
[13] in 1870. 

Higher-dimensional generalizations of CPT [d > 3) were obtained quite re¬ 
cently in and the related Cayley-type conditions were derived by the au¬ 
thors [Hj. 

The main goal of this section is to present detailed proof of Cayley-type 
condition for generalized CPT, together with discussions and examples. 

Consider an ellipsoid in 


X 


2 

1 


ai 


a;2 

H—^ = 1, ai > • • • > Od > 0, 

ad 


and the related system of Jacobian elliptic coordinates (Ai,..., Ad) ordered by 
the condition 

Ai > A2 > ■ • ■ > Ad. 


If we denote: 


Qx{x) = 


-f •••-b 


ai — A ad — \ 

then any quadric from the corresponding confocal family is given by the equation 
of the form: 


Qx : Qx{x) = 1. 


(4) 


The famous Chasles theorem states that any line in the space is tangent 
to exactly d — 1 quadrics from a given confocal family. Next lemma gives an 
important condition on these quadrics. 


Lemma 6 Suppose a line £ is tangent to quadrics Qai, ■ ■ •, from the 

family Then Jacobian coordinates (Ai,..., Ad) of any point on £ satisfy the 

inequalities 7^(As) >0, s = 1 ,..., d, where 

V{x) = {ai- x)... {ad - x){ai - x)... (od-i - x). 
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Proof. Let a; be a point of i, (Ai,...,Ad) its Jacobian coordinates, and y a 
vector parallel to i. The equation Q\{x + ty) = 1 is quadratic with respect to 
t. Its discriminant is: 


^\{x,y) = Qx{x,yf - Q\{y){Q\{x) - l), 


where 


By m. 


Q\{x,y) 


xiyi 

(X\ — A 


XdVd 

(Id A 


^\{x,y) 


(oi — A)... {ad-i — A) 
(ai - A)... (od - A) 


For each of the coordinates A = As, (l<s<(i), the quadratic equation has a 
solution t = 0; thus, the corresponding discriminants are non-negative. This is 
obviously equivalent to V{Xs) >0. □ 


4.1 Billiard inside a Domain Bounded by Confocal Qua¬ 
drics 

Suppose that a bounded domain C is given such that its boundary dfl lies 
in the union of several quadrics from the family (|1]). Then, in elliptic coordinates, 
n is given by: 

/3;<Ai</3", ..., P'd<Xd<Pl 

where Og+i < /3' < /?" < Og for 1 < s < d — 1 and —oo < < /3’f < ad- 

Consider a billiard system within fl and let Qai, • ■ •, Qad-i be caustics of 
one of its trajectories. For any s = 1,..., d, the set Ag of all values taken by 
the coordinate Ag on the trajectory is, according to Lemma 6 , included in Ag = 
{As [/Jg, /3g ] : V{X) > 0 }. By [53], each of the intervals (og+i, Og), (2 < s < d) 
contains at most two of the values Oi,..., ad-i, the interval (— oo, ad) contains 
at most one of them, while none is included in (oi, -l-oo). Thus, for each s, the 
following three cases are possible: 

First case: ai,aj G [/J',/?"], at < aj. Since any line which contains a segment 
of the trajectory touches and Qaj, the whole trajectory is placed between 
these two quadrics. The elliptic coordinate Ag has critical values at points 
where the trajectory touches one them, and remains monotonous elsewhere. 
Hence, meeting points with with Qa. and Qaj are placed alternately along the 
trajectory and Ag = Ag = [ai,aj]. 

Second case: Among ai,... ,ad-i, only ai is in [/Jg,/?"]- V is non-negative 
in exactly one of the intervals: [ag+i,ai], [ai,ag], let us take in the first one. 
Then the trajectory has bounces only on If Ui ^ P”, the billiard particle 
never reaches the boundary QfS". The coordinate Ag has critical values at meet¬ 
ing points with and the caustic Qa-, and remains monotonous elsewhere. 
Hence, As = A'^ = [/3g,ai]. If V is non-negative in [ai,ag], then we obtain 
Ag = A' =K,/3"]. 
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Third case: The segment [/3',/3"] does not contain any of values ai, ..., ad-i- 
Then V is non-negative in The coordinate As has critical values only 

at meeting points with boundary quadrics Q/j/ and 2 / 3 ", and changes monoto¬ 
nously between them. This implies that the billiard particle bounces of them 
alternately. Obviously, As = = [/3',/3"]. 

Denote [ 7 s, 7s1 := = AJ,. Notice that the trajectory meets quadrics of any 

pair Q^i , Q^n alternately. Thus, any periodic trajectory has the same number 
of intersection points with each of them. 

Let us make a few remarks on the case when 7 ' = Os+i or 7 '' = Og. This 
means that either a part of dVL is a degenerate quadric from the confocal family 
or D is not bounded, from one side at least, by a quadric of the corresponding 
type. Discussion of the case when D is bounded by a coordinate hyperplane does 
not differ from the one we have just made. On the other hand, non-existence of 
a part of the boundary means that the coordinate As will have extreme values at 
the points of intersection of the trajectory with the corresponding hyperplane. 
Since a closed trajectory intersects any hyperplane even number of times, it 
follows that the coordinate As is taking each of its extreme values even number 
of times during the period. 


Theorem 2 A trajectory of the billiard system within D with caustics Qai, ■ ■ ■, 
Qaa-i periodic with exactly Ug points at Q.y/ and Ug points at (1 < s < d) 
if and only if 

d 

{AiPyJ-A{PYj))=0 

S^l 

on the Jacobian of the curve 

r : = 'P{x) := (oi - x) • • • (a^ - x){ai - x) ■ ■ ■ {ad-i - x). 


Here, A denotes the Abel-Jacobi map, where P^i, P~f" are points on T with 
coordinates Py^ = ( 7 ', (- 1 )V^( 7 D) , -P 7 " = {is. (-l)V^( 7 s)) • 

Proof. Following Jacobi m and Darboux m, let us consider the equations: 


E 


d\s 


0. E 


XsdXg 

JW7) 


= 0, 


As ^d\g 


= 0, 


(5) 


where, for any fixed s, the square root yJV{Xg) is taken with the same sign in 
all of the expressions. Then ([5]) represents a system of differential equations of 
a line tangent to Qai, ■ ■ •, Qa^-i- Besides that. 


E 


Xt^dXg 


2d£, 


where d£ is the element of the line length. 


( 6 ) 
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Attributing all possible combinations of signs to ..., we 

can obtain 2 ‘^~^ non-equivalent systems which correspond to 2 '^ ^ different 
tangent lines to Qai, • ■ •, Qad~i from a generic point of the space. Moreover, 
the systems corresponding to a line and its reflection to a given hyper-surface 
As = const differ from each other only in signs of the roots y/V{Xs). 


Solving (O and (O as a system of linear equations with respect to 


we obtain: 


dXs 2 d£ 


dXs 




Thus, along a billiard trajectory, the differentials (—1)'* ^==== stay always 

positive, if we assume that the signs of the square roots are chosen appropriately 
on each segment. 

From these remarks and the discussion preceding this theorem, it follows that 
A* dX 

the value of the integral f — ^ between two consecutive common points of 

V'PiXs) 

the trajectory and the quadric (or Q-y'j) is equal to: 




XldX, 


'7s 


+\/V(Xs 


Now, if p is a finite polygon representing a billiard trajectory and having 
exactly points at Qj> and at (1 < s < d), then 




XldXs ,, 

/ - , . , , (1 < I < d). 


\/'P(Xs) ^ As +'\/7^(As) 

Finally, the polygonal line is closed if and only if 


which was needed. 


X^dX, ^ (1 < i < d - 1), 




□ 


Example 2 Consider two domains fl' and fl" in R^. Let 17' be bounded by the 
ellipsoid Qo and the two-folded hyperboloid Qp, 02 < /3 < ai, in such a way 
that 17' is placed between the branches of Qp. On the other hand, suppose 17" 
is bounded by Qo, the righthand branch of Qp (this one which is placed in the 
half-space xi > 0) and the plane xi = 0. Elliptic coordinates of points inside 
both 17' and 17" satisfy: 


0 ^ A3 ^ U3, ^ Ai ^ ui- 

Consider billiard trajectories within these two domains, with caustics and 
2^12 1 0,3 < Pi < 02 , a 2 < P 2 < ai- Since p 2 < /3, the segments (s S {1, 2, 3}) 
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of all possible values of elliptic coordinates along a trajectory are, for both 
domains: 

Ai = [/3,oi], A2 = [/ri,a2], A3 = [0,03]. 

In Q.", existence of a periodic trajectory with caustics and which 

becomes closed after n bounces at Qg and 2m bounces at Qp is equivalent to 
the equality: 


n(^(Po) - A{Pa,)) + 2 m{A{Pp) - A{P^,)) = 0, 

on the Jacobian of the corresponding hyperelliptic curve. In 11', existence of a 
trajectory with same properties is equivalent to: 

n(^(Po) —-^(Poa)) + 2m(^(P/3) — J) = 0 and n is even. 




Figure 3: A trajectory inside fl' Figure 4: A trajectory inside fl" 

The fact that the second equality implies the first one is due to the following 
geometrical fact: any billiard trajectory in can be transformed to a trajectory 
in Vt" applying the symmetry with respect to the a;2a:3-plane to all its points 
placed under this plane. Notice that this correspondence of trajectories is 2 to 
1 “ in such a way, a generic billiard trajectory in corresponds to exactly 
two trajectories in II'. An example of such corresponding billiard trajectories is 
shown on Figures [3] and [H 

4.2 Billiard Ordered Game 

Our next step is to introduce a notion of bounces “from outside” and “from 
inside”. More precisely, let us consider an ellipsoid Q\ from the confocal family 
0 such that A G (us+i, Qs) for some s G {!,..., c?}, where Ud+i = —00. 

Observe that along a billiard ray which reflects at Q\, the elliptic coordinate 
Ai has a local extremum at the point of reflection. 

Definition 1 A ray reflects from outside at the quadric Qx if the reflection 
point is a local maximum of the Jacobian coordinate As, and it reflects from 
inside if the reflection point is a local minimum of the coordinate Ag. 
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Figure 5: Reflection from outside Figure 6: Reflection from inside 

On Figures O and El reflections from inside and outside an ellipse and a 
hyperbola are sketched. 

Let us remark that in the case when Q\ is an ellipsoid, the notions introduced 
in Definition 1 coincide with the usual ones. 

Assume now a fc-tuple of confocal quadrics . •., from the confocal 
pencil (H]) is given, and (ii, ... ,ik) G { — 1,1}^- 

Definition 2 The billiard ordered game joined to quadrics , ■ ■ ■, 2/3^ , with 
the signature (ii,..., ik) is the billiard system with trajectories having bounces 
at , • ■ •, QfSk respectively, such that 

the reflection at is from inside if ig = +1; 
the reflection at is from outside if is = —1- 

Note that any trajectory of a billiard ordered game has d — 1 caustics from 
the same family ®. 

Suppose 2/3^, ..., 2/3fc are ellipsoids and consider a billiard ordered game 
with the signature (ii,... ,ifc). In order that trajectories of such a game stay 
bounded, the following condition has to be satisfied: 

is — 1 is+i — is—i — 1 and /Is+i ^ /Is; Ps—i ^ Ps- 

(Here, we identify indices 0 and fc + 1 with k and 1 respectively.) 

Example 3 On Figured a trajectory corresponding to the 7-tuple 

(2l, 22; 2l; 23; 22; 23; 2l); 

with the signature (1, —1,1, —1,1,1,1), is shown. 

Theorem 3 Given a billiard ordered game within k ellipsoids 2/3i; • ■ ■; 2/3^ 
with the signature (ii,...,ik)- Its trajectory with caustics Qai, ■ ■ ■, Qad,_i is 
k-periodic if and only if 

k 

Y,is{A{Pp^) - AiPof)) 
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Figure 7: Billiard ordered game 


is equal to a sum of several expressions of the form: yA{Paj,) — A{Pa^,)j on 

the Jacobian of the curve T : = V{x), where P^^ = ^/3s, +-\/7^(/3s )ct = 

minjad, ai,..., a^-i} and Qa^, Qa^, pairs of caustics of the same type. 

When Q/ 3 ^ = • • • = and ii = • • • = = 1 we obtain the Cayley-type 

condition for the billiard motion inside an ellipsoid in 

We are going to treat in more detail the case of the billiard motion between 
two ellipsoids. 


Proposition 3 The condition that there exists a closed billiard trajectory be¬ 
tween two ellipsoids and which bounces exactly m times to each of 

them, with caustics Qai, ■ • ■, Qa^-n *s- 


rank 


/ fiiPh.) 


mpj.) 


fL-d+liPh2) 

(P/3.) 


\ 


< TO — d -I- 1. 






Here 


fj 


y- Bq- Bi{x -Pi) -Pd+j- 2 (x - ^ 

xd+j-l 


l<j<m — d-\-l, 


and y = Bq-\-Bi{x—(3i)-\-. .. is the Taylor expansion around the point symmetric 
to Pfj.^ with respect to the hyperelliptic involution of the curve F. (All notations 
are as in Theorem [3l) 


Example 4 Consider a billiard motion in the three-dimensional space, with 
ellipsoids Qo and as boundaries (0 < 7 < 03 ) and caustics Qai and Qa.- 
Such a motion closes after 4 bounces from inside to Qo and 4 bounces from 
outside to Q.y if and only if: 

rankAl < 2 . 
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The matrix X is given by: 

Xn = -SCo + C^-f + 3Bo + 2 Bi7 + 

Xi2 = —4Co + Ci7 + 4i?o + 3 -Bi7 + 2B2"f^ + B^j^ 

X 21 = 6 C 0 - 3(717 + C27 " - 6^0 - 3Si7 - B 27" 

X22 = 10(7o - 4(7 i 7 - lOBo - 6S17 - 3^27^ - ^37^ 

X31 = -10(7o + 6(7 i 7 - 3(727^ + <7372 + IOBq + 4Si7 + ^27^ 

X32 = -20(7o - 10(7 i 7 - 4(727^ + <^37^ + 20Bo + IOS17 + 4^27^ + S37^ 
and the expressions 

— \/(oi — x){a 2 — x){a 3 — x){ai — x)(a 2 — x) = i?o + Bix + B 2 X^ + ..., 

+ \/(oi — x){a 2 — x){a 3 — x){ai — x){a 2 — x) = Cq + Ci{x — ^) + ... 

are Taylor expansions around points x = 0 and x = j respectively. 

Example 5 Using the same notations as in the previous example, let us con¬ 
sider trajectories with 4 bounces from inside to each of Qo and as shown 
on Figure m 



Figure 8: A closed trajectory of the billiard 
ordered game with 8 alternate bounces from 
inside of two ellipses 


The explicit condition for periodicity of such a trajectory is: 

rankA < 2, 


with 

All = -4(7o -b (7 i 7 -b 3^17 -b 2^27^ -b B^j^ 

Ai 2 = -3(7o -b (7 i 7 -b 3Bo + 2 Bij + ^27^ 

Ai2i = -6(7o -b (727^ -b 6B0 + 6 B 1 J + 5^27^ + 3^37^ 
A22 = -6(7o -b (7 i 7 -b (727^ + 6S0 -b 6^17 -b 3^27^ 
A31 = —4(7o -b (737^ -b 450 -b 45 i 7 -b 4527^ -b 3^37^ 
A32 = -4(7o -b (727^ + (737^ + 45o -b 45i7 -b 3527^. 
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5 Elliptical Billiard as a Discrete Time System: 
d — Arbitrary, k = 1 

Another approach to the description of periodic billiard trajectories is based on 
the technique of discrete Lax representation. 

In this section, first, we are going to list the main steps of algebro-geometric 
integration of the elliptic billiard, following [32] ■ Then, the connection between 
periodic billiard trajectories and points of finite order on the corresponding 
hyperelliptic curve will be established, and, using results from Section [31 the 
Cayley-type conditions will be derived, as they were obtained by the authors 
in |191120| . In addition, here we provide a more detailed discussion concerning 
trajectories with the period not greater than d and the cases of the singular 
isospectral curve. 

Following [321, the billiard system will be considered as a system with the 
discrete time. Using its integration procedure, the connection between periodic 
billiard trajectories and points of finite order on the corresponding hyperelliptic 
curve will be established. 


5.1 XYZ Model and Isospectral Curves 

Elliptical Billiard as a Mechanical System with the Discrete Time. 

Let the ellipsoid in be given by 

{Ax, x) = 1. 

We can assume that A is a diagonal matrix, with different eigenvalues. The 
billiard motion within the ellipsoid is determined by the following equations: 


— l^kyk-\-l 


where 


yfc + l Uk — 


i^k') 

{Ayk+i,yk+iy 


‘2{Axk , yk) 
{Axk,Axk)' 


Here, Xk is a sequence of points of billiard bounces, while yk 
are the momenta. 


^k ^k—1 
\^k ^k—1 1 


Connection between Billiard and XYZ Model. To the billiard system with 
the discrete time, Heisenberg XYZ model can be joined, in the way described 
by Veselov and Moser in [32| and which is going to be presented here. 

Consider the mapping ip : {x, y) {x', y') given by 

x'f, = Jyk+i= J{yk+VkAxk), y'^, = -J~^Xk, J = A~^. 
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Notice that the dynamics of tp contains the billiard dynamics: 

x'k = Jy'k+i = -Xk+i-, y'k = = -yk+i, 

and define the sequence {xk,yk)' 

{xo,yo) := {xo,yo), {xk+i,yk+i) ■= (p{xk,yk), 

which obeys the following relations: 

Xk-\-l — Jyk ^kJ Xk-i yk+1 — J Xkt 

where the parameter i/k is such that \yk\ = 1, {Axk,Xx) = 1- This can be 
rewritten in the following way: 

Xk-\-l “f Xk — 1 — k'kJ Xk- 

Now, for the sequence qk ■= J~^Xk, we have: 

Qk+i + Qk-i =’^kJ~^qk, |gfc| = l- 

These equations represent the equations of the discrete Heisenberg XYZ system. 

Theorem 4 [32j Let {xk,yk) be the sequenee eonnected with elliptical billiard 
in the described way. Then qk = J~^Xk is a solution of the discrete Heisenberg 
system. 

Conversely, if qk is a solution to the Heisenberg system, then the sequence 
Xk = (—l)^J 92 fc is a trajectory of the discrete billiard within an ellipsoid. 

Integration of the Discrete Heisenberg XYZ System. Usual scheme 
of algebro-geometric integration contains the following [32]. First, the sequence 
Lfe(A) of matrix polynomials has to be determined, together with a factorization 

L(A) = B{X)C{X) ^ C{X)B{X) = B'{X)C'{X) = L'{X), 

such that the dynamics L L' corresponds to the dynamics of the system qk. 
For each problem, finding this sequence of matrices requires a separate search 
and a mathematician with the excellent intuition. All matrices Lk are mutually 
similar, and they determine the same isospectral curve 

F : det(L(A) — pi) = 0. 

The factorization Lk = BkCk gives splitting of spectrum of Lk. Denote by tpk 
the corresponding eigenvectors. Consider these vectors as meromorphic func¬ 
tions on F and denote their pole divisors by Dk. 

The sequence of divisors is linear on the Jacobian of the isospectral curve, 
and this enables us to find, conversely, eigenfunctions ipk, then matrices Lk, 
and, finally, the sequence (qk)- 
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Now, integration of the discrete XYZ system by this method will be shortly 
presented. Details of the procedure can be found in [3^ . 

The equations of discrete XYZ model are equivalent to the isospectral de¬ 
formation: 

Lfe+i(A) =Afe(A)Lfe(A)A-i(A), 

where 

LkW = + Mk-l JQk — A^q/c-l 0 Qk-lt 

AkiX) = J - XqkZ) Qk-i- 

The equation of the isospectral curve T : det(L(A) — fil) = 0 can be written in 
the following form: 

d—1 d 

i=i j=i 

where v = XYli^iil^ — l^i) and /ri ,..., fJ-d-i are zeroes of the function: 




E 


Fi{x,y) 


F,=xf 


{x AJy) 


3=Ai 


Jf - Jf 


X = qk-1, y = Qk- 


It can be proved that y,i,, jid-i are parameters of the caustics corresponding 
to the billiard trajectory |33j . Another way for obtaining the same conclusion is 
to calculate them directly by taking the first segment of the billiard trajectory 
to be parallel to a coordinate axe. 

If eigenvectors il^k of matrices Lk{X) are known, it is possible to determine 
uniquely members of the sequence qk. Let Dk be the divisor of poles of function 
4’k on curve T. Then [5^ : 


Dk+i — F>k + Poo — Pq, 


where Poo is the point corresponding to the value /i = oo and Pq to y, = 0, 

A — (.Qk'! J Qk+d 


5.2 Characterization of Periodical Billiard Trajectories 

In the next lemmae, we establish a connection between periodic billiard se¬ 
quences Qk and periodic divisors Dk. 

Lemma 7 [19] Sequence of divisors Dk is n-periodic if and only if the sequence 
qk is also periodic with the period n or Qk+n = ~Qk for all k. 
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Proof. If, for all k, qk+n = qk, or qu+n = —qk for all /c, then, obviously, 
Lk+n = Lk- Thus, the sequence of eigenvectors tpk is periodic with the period 
n. It follows that the sequence of divisors is also periodic. 

Suppose now that Dk+n = Dk for all k. This implies that ifk+n = Ckipk- We 
have: 

= Ak{X)i>k. 

Let fj,i and fj ,2 be values of parameter ^ which correspond to the value A = 1 on 
the curve T, and 

'Sk = ('0fc(l,Mi),V'fe(l,M2))- 

From V'fc+i = Ak{X)'4>k , we obtain Ak{l) = 'I'fc+i'I'^^. It follows that 

Ak{l) = 

From the condition det Afe = det Afe+i for all fc, we have Ck = Ck+i- Thus, the 
sequence 

Afe(l) = J - qk® qk-i 

is n-periodic. From there. 


1 

Q/c+n — OikQk-, Q_k-\-n—l — Q_k—1- 

Otk 


Since |gfc| = 1, we have offc = 1 or afc = — 1, where all ak are equal to each other, 
which proves the assertion. □ 


Lemma 8 [12] The billiard is, up to the central symmetry, periodic with the pe¬ 
riod n if and only if the divisor sequence Dk joined to the corresponding Heisen¬ 
berg XYZ system is also periodic, with the period 2n. 

Proof. Let Xk+n = ctXk for all k, a G { — 1,1}. Join to a billiard trajectory 
{xk,yk) the corresponding flow {xk,yk). Since 

{x2k,y2k) = {-l)'"{xk,yk), {x2k+i,y2k+i) = (I}{x2k,y2k), 


where the mapping (j) is linear, we obtain: 

^k-\-2n — 0 :( 1 ) Xk. 

From there, immediately follows that qk+ 2 n = (^{—Yj'^qk. According to 
Lemma 3, the divisor sequence Dk is 2n-periodic. □ 

Applying the previous lemma, we obtain the main statement of this section: 

Theorem 5 [20] A condition on a billiard trajectory inside ellipsoid Qq in 
with non-degenerate caustics Q^iii ■ ■ ■, Q^id-n periodic, up to the central 

symmetry, with the period n > d is: 



/ Bn+l 

Bn 

■ Bd+1 \ 

rank 

Bn-\-2 

Bn-\-l 

Bd+2 


\ B2n-1 

B2n-2 
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where 


y/{x- ^li)---{x- lJ.d-l){x - Oi) • • • (x - Qd) = So + BiX + B 2 x‘^ + - 

Proof. The trajectory is periodic with period n if, by LemmalU the correspond¬ 
ing divisor sequence on the curve T has the period 2n, i.e. 2n(Poo ~ Bo) = 0 on 
J'(r). Curve r is hyperelliptic with genus g = d — 1. Taking A{Poo) to be the 
neutral on J"(r) we get the desired result by applying Lemma O □ 

Cases of Singular Isospectral Curve. When all Oi,..., a^,/ii, ■ ■ ■, ^d-i 
are mutually different, then the isospectral curve has no singularities in the 
affine part. However, singularities appear in the following three cases and their 
combinations: 

(i) Qi = fij for some i,j- The isospectral curve (O decomposes into a rational 
and a hyperelliptic curve. Geometrically, this means that the caustic corre¬ 
sponding to g,i degenerates into the hyperplane Xi = 0. The billiard trajectory 
can be asymptotically tending to that hyperplane (and therefore cannot be pe¬ 
riodic), or completely placed in this hyperplane. Therefore, closed trajectories 
appear when they are placed in a coordinate hyperplane. Such a motion can be 
discussed like in the case of dimension c? — 1. 

(ii) Gi = Gj for some i ^ j. The boundary Qq is symmetric. 

(iii) g-i = gj for some i ^ j. The billiard trajectory is placed on the corre¬ 
sponding confocal quadric hyper-surface. 

In the cases (ii) and (iii) the isospectral curve T is a hyperelliptic curve with 
singularities. In spite of their different geometrical nature, they both need the 
same analysis of the condition 2nPo ~ 2nE for the singular curve ©. 

Immediate consequence of Lemma |4] is that Theorem [5] can be applied not 
only for the case of the regular isospectral curve, but in the cases (ii) and (iii), 
too. Therefore, the following interesting property holds. 

Theoreui 6 If the billiard trajectory within an ellipsoid in d-dimensional Eu- 
cledean space is periodic, up to the central symmetry, with the period n < d, then 
it is placed in one of the n-dimensional planes of symmetry of the ellipsoid. 

Proof. This follows immediately from Theorem [5] and the fact that the section 
of a confocal family of quadrics with a coordinate hyperplane is again a confocal 
family. □ 

Note that all trajectories periodic with period n up to the central symmetry, 
are closed after 2n bounces. A statement sharper than Theorem [6] for the tra¬ 
jectories closed after n bounces, where n < d, can be obtained in the elementary 
fashion: 

Proposition 4 If the billiard trajectory within an ellipsoid in d-dimensional 
Eucledean space is periodic with period n < d, then it is placed in one of the 
(n — 1)-dimensional planes of symmetry of the ellipsoid. 
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Proof. First, consider the case n = d. Let xi ■ • ■ Xdhe a periodic trajectory, and 
{N, x) = a the equation of the hyperplane spanned by its vertices. Here, is a 
vector normal to the hyperplane and a is a constant. Since all lines normal to 
the surface of the ellipsoid at the points of reflection belong to this hyperplane, 
it follows that {AN,Xi) = {N,Axi) = 0. Thus, {AN,x) = 0 is also an equation 
of the hyperplane, so a = 0 and the vectors N, AN are collinear. From here, 
the claim follows immediately. 

The case n < d can be proved similarly, or applying Theorem [5] and the 
previous case of this proposition. □ 

This property can be seen easily for d = 3. 

Example 6 Consider the billiard motion in an ellipsoid in the 3-dimensional 
space, with /ri = ^ 2 , when the segments of the trajectory are placed on gen¬ 
eratrices of the corresponding one-folded hyperboloid, confocal to the ellipsoid. 
If there existed a periodic trajectory with period n = d = 3, the three bounces 
would have been coplanar, and the intersection of that plane and the quadric 
would have consisted of three lines, which is impossible. It is obvious that any 
periodic trajectory with period n = 2 is placed along one of the axes of the 
ellipsoid. So, there is no periodic trajectories contained in a confocal quadric 
surface, with period less or equal to 3. 


6 Periodic Trajectories of Billiards on Quadrics 


in Rd 

In [12], the billiard systems on a quadric Qo in 



0± ^ ^ CLdj 


are defined as limits of corresponding billiards within Qo, when one of the 
caustics tends to Qq. The boundary of such a billiard consists of the intersection 
of Qo with certain confocal quadrics Q/?,,..., Q/j^ of the family ([4|). Between 
the bounces, the billiard trajectories follow geodesics of the surface of Qo, while 
obeying the reflection law at the points of the boundary. A tangent line issued 
from any point of a given trajectory touches, besides Qo, also d — 2 confocal 
quadrics Qai, ■ ■ ■, from the family (l4|). These d — 2 quadrics are fixed 

for each trajectory, and we shall refer to them as caustics. 

The question of description of periodic trajectories of these systems was 
formulated as an open problem by Abenda and Fedorov [1]. 

Like in Section [d] we will consider two cases: the billiard inside a domain 
H C Qo bounded by several confocal quadrics, and the billiard ordered game 
within several confocal quadrics of the same type. 

For the first case, the domain 17 is given by: 


/ 3 - 


[// 


, Pd-l — ^d-l ^ Pd-li iPd — 0), 
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where /Jg,/?" G [as+i,as] for 1 < s < d — 1. Like in Section [H Ag, (1 < s < 
d — 1), is defined as the set of all values taken by the coordinate Ag on the given 
trajectory, and it can be shown in the same way that 

Ag = {A G [/?:,/?"] : 7^i(A)>0}, 

with 

Viix) = -x{ai - x) ■■■ {ad - x){ai - x) ■ ■ ■ {ad -2 - x), 

where oi, ..., ad -2 are parameters of caustics of the trajectory. 

Denote [ 7 ', 7 "] — A*. 

Before formulating the theorem, let us define the following projection of the 
Abel-Jacobi map on the curve 




/ 

Theorem 7 A trajectory of the billiard system constrained on the ellipsoid Qo 
within fl, with caustics Qai, ■ ■ ■, Qai_ 2 ! periodic with exactly Ug bounces at 
each of quadrics Q^'^, Q-y'j, (l<s<d—1), if and only if 


Li : y^=Vi{x), 


by: 


A{P) = 


^ xdx 

Jo — 

fP x^dx 

Jo 


f. 


.p X 


y 


d-l 

J2ns{A{PYj-A{PYj))=0. 

S^l 

Here P-y', P^ij are the points onTi with coordinates P^/^ = ^ 7 ', (—l)®-\/Pi( 7 g)^ , 

PYj = (isA-^rV'p^)- 

Proof. By [27] , the system of differential equations of a geodesic line on Qo with 
the caustics Q^, ..., Qad ,_2 is: 


d-l 


E 


XgdXs 


0^ E 


A^dAg 


= 0 , 


Y Af-"dAg 


= 0 , 


with the square root \/P{Xs) taken with the same sign in all of the expressions, 
for any fixed s. Also, 


^ Af-^dAg 


= 2de, 
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where M is the length element. 

Now, the rest of the demonstration is completely parallel to the proof of 
Theorem [5] □ 

In the same way as in Section 31 a billiard ordered game constrained on the 
ellipsoid Qo within given quadrics , ■ • ■, of the same type can be defined. 
The only difference is that now the signature a = (ii,...,ife) can be given 
arbitrarily, since trajectories are bounded, lying on the compact hypersurface 
Qo- Denote by Qai, ■ ■ ■, caustics of a given trajectory of the game. 

Since quadrics ,. ■ •, are all of the same type, there exist /i', p." in the set 
S = {oi,..., Od, ai,..., ad- 2 } such that /3i,...,/3fe S [^', p,"] and (/j,', n S' is 
empty. 

Associate to the game the following divisors on the curve T 1 : 

if is = Zs-i-i = 1 
if is = -is+i = l,Ps < Ps+i 
or is = -is+i = -1, Ps > Ps+i 
if is — is+i — I 5 Ps Ps+i 
if Is = -is+i = -1, A < Ps+i 
if fs — fs+l — 1 ; 

where and P^" are the branching points with coordinates (/i', 0 ) and (/r", 0 ) 
respectively. 

Theorem 8 Given a billiard ordered game constrained on Qo within quadrics 
Q/ 3 i; • • ■: Q/ 3 fc with signature a = (*i,... , 1 ^). Its trajectory with caustics Qai, 

■ ■ ■! Qad -2 k-periodic if and only if 

k 

'£^s{A{P0P - A{Vs)) 

S = 1 

is equal to a sum of several expressions of the form A^PaA — A{Pa^,) on the 
Jacobian of the curve Ti y^ = Vi{x), where Pp^ = ^Ps, +a/Pi(Ps)^ and Qa^, 
Qdp, ore pairs of caustics of the same type. 

Example 7 Consider the case d = 3 and a billiard system constrained on the 
ellipsoid Qo with the boundary and caustic Qc, 03 < 7 < a < 02 . A 
sufficient condition for a corresponding billiard trajectory to be n-periodic is: 

n(A(P.,)-A(P„)) =0, 

or, equivalently, in Cayley-type form: 



/ 

Cp+i 

Cp+2 

1 0 

to 

1 

to 




rank 


Cp+2 

Cp+3 

(^2p-l 


<p-2, 

cu 

II 

e 


1 

C2p 

C 2 P+I 

Csp-s 





Vs = 


fl" 5 


P^l" - P/i', 
Pfj.' - Pfi ", 

Pc, 
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and 



/ 

Cp+1 

Cp+2 

C2P-I 

rank 


Cp+2 

Cp+3 

■ C2p 


V 

C2p 

C2P+I 

C3P-2 


n = 2p + 1, 


where 


2/ = Co + Cl 



+ C2 




is the Taylor expansion around the point P™, with x = -. (See also 1211.1 

a — X 

More precisely, this condition will be satisfied if and only if the trajectory is 
n-periodic and its length L with respect to the parameter s: 


ds = X 1 X 2 ■■ ■ Xndl 


is such that the vector 


/ L/2 \ 

0 


V 0 


+ nA{Pj) — nA{Pa) 


belongs to the period-lattice of the Jacobian of the corresponding hyper-elliptic 
curve. 


Example 8 In a recent preprint [T] , another interesting class of periodical bil¬ 
liard trajectories on 2-dimensional ellipsoid, associated with closed geodesics, is 
described. Such geodesics are closely connected with 3-elliptic KdV solutions, 
which historically goes back to Hermite and Halphen |26j . 

The corresponding genus 2 isospectral curve is: 

G : y'^ = -^{‘ix^-9g2X-27g3){x'^-ig2). (8) 

This curve is 3 : 1 tangetial covering of two elliptic curves: 

El : yl= 4xf - g 2 Xi - gs, 

E 2 : yl = xl-‘^{gl + 9gl)x2-‘^i3gl-gsgl). 

The corresponding conditions for closed billiard trajectories are obtained in [1], 
by use of Lemma 3 from |20j , which corresponds to Lemma [3] of the present 
paper. For the exact formulae and details of calculations, see [I]. 


Note that in the same algebro-geometric situation as in the previous example, 
with the application of Theorems [7] and [3 it is possible to describe a wider class 
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of closed billiard trajectories, without the assumption that their segments belong 
to closed geodesics. Namely, an important property of the covering G —> Ei from 
Example [5] is that the holomorphic differential on G reduces to — | on 
El. Thus, whenever the isospectral curve can be represented in the form ([8]), it 
is possible to obtain conditions of Gayley’s type by the application of Theorem 
[71 Theorem |S] and Lemma [31 

This is going to be illustrated in the next example. 

Example 9 Consider the billiard motion on quadric Qq in with oi = oa/S, 
02 = |a, 03 = —U'/S, with the caustic Qa, a = — |o, o > 0. The domain fi is 
bounded by confocal ellipsoid Qp, j3 < By Theorem [71 the condition for 

such a billiard to be n-periodic is: 

n{A{P0) - A{P_^^)) = 0 (9) 

on the Jacobian of the curve 

= —X ~ 3a^). 

This curve is of the form [51 with ^2 = ffs = 0. 

The covering G ^ Ei is given by: 

1 2 x^ — 9o^x 

^“9x2-3o2’ ^ ^ (x2 - 3o2)2 ' 

The condition @ reduces to the following relation: 

n[A{Qp) - ^(Qoo)) = 0 

on the Jacobian of the curve: 


El ■. yl = 4xf - o^xi. 


with the point Qp G Ei having the coordinates {zi = (}, Wi = 4/3^ — a^/3), 
1 /3^ 

/3 = — 9 Q<x> being the infinite point on Ei. 

9 — 3a2 

By Lemma |3l the desired condition becomes: 


^3 

Ba 

Bjn+l 




Ba 

B5 

Bm+2 

= 0, 

for n 

= 2m 

-^m+1 

Bm+2 

B2m-1 




B2 

B3 

Bm+1 




Bs 

Ba 

Bm+2 

= 0, 

for n 

= 2m + 1 

Bm-\-l 

Bm+2 

B2m 





where \/4xf — a2xi = i3o+J3i(xi—/3)+i32(xi—/3)^ + ... is the Taylor expansion 
around the point Qp. 
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We will finish this section by analysis of the behaviour of geodesic lines on 
quadrics after the reflection of a confocal quadric. 

Denote by g a geodesic line on the quadric Qo in and by D a domain 
on Qo bounded by a single confocal quadric Qp. We will suppose that the set 
g n Q^ ^ 0 and that the geodesic g intersects Q^ transversally. Under these 
assumptions, the number of their intersection points will be finite if and only if 
the line g is closed. These points are naturally divided into two sets - one set, 
denote it by Si, contains those points where g enters into U, while S 2 contains 
the points where the geodesic leaves the domain. 

Consider reflections on Qjj according to the billiard rule in each point of 82 - 
Applying these reflections to g, we obtain a family of geodesic segments on Qq. 
It appears that all these segments belong to one single geodesic line, see Figure 

H 

Proposition 5 Let g' be a geodesic line which contains a point P G S 2 , such 
that g and g' satisfy the reflection law in P at Qp. Then g' contains all points 
of S 2 and the two lines g, g' satisfy the reflection law at Q^g in each of these 
points. 

If g is closed, then g' is also closed and they have the same number of 
intersection points with Q^. 



Figure 9: Reflection of a geodesic line 


Proof. The reflection from inside at Q^g corresponds to the shift by the vector 
A{Pp) — A{tP 0 ) on the Jacobian of the curve Fi, where t : Fi ^ Fi is the 
hyper-elliptic involution: {x,y) 1 -^ (x,—y). Thus, the shift does not depend on 
the choice of the intersection point and it maps one geodesic into another one. 
This proof is essentially the same as the proof of Lemma 5.1 in [1], which covers 
the special case of a closed geodesic line on the two-dimensional ellipsoid. □ 

Observe that g' not necessarily contains the points of Si, see Figure [9] 

Denote by Tfs the above mapping g 1 —> g' in the set of geodesics on Qo (on 
geodesic lines that are tangent to Qp or do not intersect it at all, we define Tp 
to be the identity). As a consequence of Proposition [5l we obtain a new porism 
of Poncelet type. 
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Corollary 2 Suppose that T^{g) = g. Then is the identity on the class of 
all geodesics sharing the same caustics with g. 

Moreover, if the geodesic line g is closed, then any billiard trajectory in O 
with the initial segment lying on g, will be periodic. 


7 Virtual Billiard Trajectories 

Apart from reflections from inside and outside of a quadric, as we defined in 
Section 4, which correspond to the real motion of the billiard particle, it is 
of interest to consider virtual reflections. These reflections were considered by 
Darboux in m The aim of this section is to prove and generalize a property of 
virtual reflections formulated by Darboux in the three-dimensional case in the 
footnote [TS] on p. 320-321: 

.. II importe de remarquer: le theoreme donne dans le texte suppose essen- 
tiellement que les cotes du polygone soient formes par les parties reelles et non 
virtuelles du rayon refiechi. II existe des polygones fermes d’une tout autre 
nature. Etant donnes, par example, deux ellipsoides homofocaux [Ef), {Ei), 
si, par une droite quelconque, on leur mene des plans tangents, on aura quatre 
points de contact oq, bg sur (Eg), ai, bi sur (Ei). Le quadrilatere agaibgbi 
sera tel que les bissectrices des angles ai, bi soient les normales de {Ei), et 
les bissectrices des angles ag, bg les normales de {Eg), mais il ne constituera 
pas une route reelle pour un rayon lumineux; deux de ses cotes seront formes 
par les parties virtuelles des rayons reflechis. De tels polygones meriteraient 
aussi d’etre etudies, leur theorie offre les rapports les plus etroits avec celle de 
I’addition des fonctions hyperelliptiques et de certaines surfaces algebriques ... ” 

More formally, we can define the virtual reflection at the quadric Q as a map 
of a ray £ with the endpoint Pg {Pg G Q) to the ray complementary to the one 
obtained from £ by the real reflection from Q at the point Pg. 

Let us remark that, in the case of real reflections, exactly one elliptic coordi¬ 
nate, the one corresponding to the quadric Q, has a local extreme value at the 
point of reflection. On the other hand, on a virtual reflected ray, this coordinate 
is the only one not having a local extreme value. In the 2-dimensional case, the 
virtual reflection can easily be described as the real reflection from the other 
confocal conic passing through the point Pg. In higher dimensional cases, the 
virtual reflection can be regarded as the real reflection of the line normal to Q 
at Po- 

From now on, we consider the n-dimensional case. 

Let points Xi,X 2 ; Yi,Y 2 belong to quadrics Qi, Q 2 of a given confocal 
system. 

Definition 3 We will say that the quadruple of points Xi, X 2 ,Yi,Y 2 consti¬ 
tutes a virtual reflection configuration if pairs of lines XiYi, X 1 Y 2 ; X 2 Y 1 , X 2 Y 2 ', 
XiYi, X 2 Y 1 ; X 1 Y 2 , X 2 Y 2 satisfy the reflection law at points Xi, X 2 of Qi and 
Yi, Y 2 of Q 2 respectively. 
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Figure 10: Virtual reflection configuration 


Now, the Darboux’s statement can be generalized and proved as follows: 

Theorem 9 Let Qai, Q \2 ^6 confocal quadrics, Xi, X 2 points on Qai and Yi, 
Y 2 on Qa 2 • If the tangent hyperplanes at these points to the quadrics belong to 
a pencil, then Xi, X 2 ,Yi,Y 2 constitute a virtual reflection configuration. 

Furthermore, if Qai and Qa 2 are ellipsoids and Ai > X 2 , then the sides of 
the quadrilateral X 1 Y 1 X 2 Y 2 obey the real reflection from and the virtual 
reflection from Qai • 

Proof. Denote by ^ 1 , ^2 and rji, 772 the tangent hyperplanes to Qai at Xi, X 2 
and to Qa 2 at Yi, Y 2 , respectively. All these hyperplanes belong to a pencil, 
thus their poles with respect to any quadric will be collinear - particularly, the 
pole P of lies on line yiV 2 - If <5 = Y 1 Y 2 fl ^ 1 , then pairs P, Q and Yi, Y 2 are 
harmonically conjugate. It follows that the lines XiYi, X 1 Y 2 obey the reflection 
law from ^ 1 . We prove similarly that all other adjacent sides of the quadrilateral 
X 1 Y 1 X 2 Y 2 obey the law of reflection on the corresponding quadrics. 

Now, suppose that Qai, Qa 2 are ellipsoids, and Ai > A 2 . This means that 
Qai is placed inside Qa 2 > thus the whole quadrilateral A 1 V 1 A 2 V 2 is inside Qa 2 - 
This means that the reflections in points Yi, Y 2 are real reflections from inside 
on Qa 2 ■ Besides, the ellipsoid Qai is completely placed inside the dihedron with 
the sides rji, rj 2 - This ellipsoid is also inside the dihedron Z(^i,^ 2 )- Since planes 
r]i i 772 are outside Z(^i,^ 2 ), it follows that points Yi, Y 2 are also outside this 
dihedron. Thus, points Yi, I 2 are placed at the different sides of each of the 
planes ^ 1 , ^ 2 , and reflections of Qai are virtual. □ 

We are going to conclude this section with the statement converse to the 
previous theorem. 

Proposition 6 Let pairs of points Xi, X 2 and Yi, Y 2 belong to confocal ellip¬ 
soids Qi and Q 2 , and let ai, a 2 , fli, P 2 be the corresponding tangent planes. If 
a quadruple Xi, X 2 ,Yi,Y 2 is a virtual reflection configuration, then planes ai, 
0 : 2 , fli, (32 belong to a pencil. 

Proof. Consider the pencil determined by ai and fli. Let af be planes of 
this pencil, tangent respectively to Qi, Q 2 at points X 2 , Yf, and distinct from 
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Qfi and (3i. By Theorem [HI the quadruple Xi,X 2 ,Yi,Y 2 is a virtual reflection 
configuration. Moreover, if denote by Ai, A 2 parameters of Qi, Q 2 and assume 
Ai > A 2 , then the sides of the quadrangle obey the reflection law at points 
^1,1^2 the virtual reflection at Since the ray obtained from XiYi 

by the virtual reflection of Qi at Xi, has only one intersection with we have 
I 2 = ^ 2 '. Points X 2 and X 2 coincide, being the intersection of rays obtained 
from YiATi and Y 2 X 1 , by the reflection at the quadric Q 2 - Now, the four tangent 
planes are all in one pencil. □ 

8 On Generalization of Lebesgue’s Proof of Cay¬ 
ley’s Condition 

In this section, the analysis of possibility of generalization of the inspiring Le¬ 
besgue’s procedure to higher dimensional cases will be given. 

A higher dimensional analogue of the crucial lemma from [31j . which is 
Lemma [I] of this article, is the following: 

Lemma 9 Let Qi, Q 2 be quadrics of a confocal system and let lines £ 1 , £2 
satisfy the reflection law at point Xi of Qi and £ 2 , £3 at Y 2 of Q 2 - Then line £\ 
meets Q 2 at point Yi and £3 meets Qi at point X 2 such that pairs of lines £ 1 , 
Y 1 AI 2 and Y 1 X 2/3 satisfy the reflection law at points Yi, X 2 of quadrics Q 2 , Qi 
respectively. Moreover, tangent planes at Xi, X 2 , Yi, Y 2 of these two quadrics 
are in the same pencil. 

This statement can be proved by the direct application of Theorem on 
virtual reflections. Nevertheless, there is no complete analogy between Lemma 
Oand the corresponding assertion in plane. Lines £ 1 , £3 and £ 2 , Y 1 Y 2 are generi- 
cally skew. Hence we do not have the third pair of planes tangent to the quadric, 
containing intersection points of these two pairs of lines. 

Nevertheless, a complete generalization of the Basic Lemma, can be formu¬ 
lated as follows: 

Theorem 10 Let X be a dual pencil of quadrics in the three-dimensional space. 

For a given quadric Tq S IF, there exist quadruples a, (3, 7, <5 of planes 
tangent to Tq, and quadrics ri,r 2 ,r 3 G £F touching the pairs of intersecting 
lines af) and "fS, aj and (35, a5 and (3'y respectively, with the tangent planes 
to Ti, r 2 , La at points of tangency with the lines, all being in one pencil A. 
Moreover, the six intersecting lines are in one bundle. 

Every such a configuration of planes a, (3, 7, 5 and quadrics Lq, Ti, r2, Ts 
is determined by two of the intersecting lines, and the tangent planes to Ti, r 2 , 
orTs corresponding to these two lines. 

Proof. Suppose first that two adjacent lines a(3 and 07 are given. There exist 
two quadrics in T touching af3 — let Ti be the one tangent to the given plane 
/i D a[3. r2 denotes the quadric that touches 07 and the given plane tt D 07. 
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The pencil A is determined by its planes /r and tt. All three lines a/3, ay, A are 
in one bundle B, see Figure ITT] 



Figure 11: Theorem fTOl 

Note that lines a/3 and ay determine plane a and that a touches a unique 
quadric Fq from T. Thus, /3 and y are determined as tangent planes to Fq, 
containing lines a/3 and ay respectively and being different from a. 

Let V be the plane from A, other than /i, that is touching Fi. We are going 
to prove that the point of tangency v with Fi belongs to y. 

Denote by <i> the dual pencil determined by quadric Fi and the degenerate 
quadric which consists of lines a/3 and y^. Since y^ G 6, these two lines are 
coplanar. 

Dual pencils T and $ determine the same involution of the pencil ay, be¬ 
cause they both determine the pair a,y and the quadric Fi is the common 
for both pencils. Since quadric Vi ^ T determines the pair tt, tt of coinciding 
planes, a quadric Q determining the same pair has to exist in $. This quadric, 
as well as all other quadrics in $, touches v and /r. Since /r, tt all belong to 
the pencil A, Q is degenerate and contains A. Since any quadric of touches 
/i, TT at points of lines a/3, ay respectively, the other component of Q also has 
to be A, thus Q is the double A. 

It follows that all quadrics of d), and particularly Fi, are tangent to at a 
point of yi^. 

Similarly, if k 7 ^ tt is the other plane in A tangent to r 2 , the touching point 
belongs to /3 and to S, the plane, other than y, tangent to Fg and containing 
the line yz^. 

Now, let us note that and <i> determine the same involution on pencil aS, 
because they both determine pair a,S and Fi belongs to both pencils. Thus, 
the common plane p of pencils A and ai5 is tangent to a quadric of iF at a point 
of aS. Denote this quadric by Fg. Similarly as before, we can prove that Fg 
is touching the line /3y and the corresponding tangent plane a is common to 
pencils A and /3y. 

Now, suppose that two non-adjacent lines a/3, yJ, both tangent to a quadric 
Fi G iF, with the corresponding tangent planes p, v, are given. Similarly as 
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above, we can prove that the plane p is tangent to a quadric from T at a point 
of a5. So, we can consider the configuration as determined by adjacent lines 
a/3, aS and planes /i, p. In this way, we reduced it to the previous case. □ 

For the conclusion, we will generalize the notion of Cayley’s cubic. 

Definition 4 The generalized Cayley curve is the variety of hyperplanes tan¬ 
gent to quadrics of a given confocal family in at the points of a given line 
L 

This curve is naturally embedded in the dual space *. 

On Figure [T^ we see the planes which correspond to one point of the line ^ 
in the 3-dimensional space. 



Figure 12: Three points of the generalized Cay¬ 
ley curve in dimension 3 


Proposition 7 The generalized Cayley curve in for d > 3 is a hyperelliptic 
curve of genus g = d — 1. Its natural realization in * is of degree 2d — 1. 

Proof. Let us consider the projection from the generalized Cayley’s curve to the 
line of the parameters of the confocal family. Since £ intersects each quadric 
twice, this is a two-folded covering, with branching points corresponding to 
quadrics touching £ and the degenerate ones. Since there is (d — 1) -I- (d -I- 1) 
such quadrics, we obtain the genus directly from Riemann-Hurwitz formula. 

Its degree is equal to the number of intersection points with a hyperplane 
in *. Such a hyperplane is a bundle of hyperplanes containing one point in 
C'^. Take P S £ to be this point. Since there are d quadrics from the confocal 
family containing P and d — 1 tangent to £, the assertion follows. □ 

Let us note that this curve is isomorphic to the Veselov-Moser isospectral 
curve ©■ Also, in the 3-dimensional case, it is isomorphic to the Jacobi hy¬ 
perelliptic curve, which was used by Darboux considering the generalization of 
Poncelet theorem. 

Further development of these ideas will be presented in the separate publi¬ 
cation [22] . 
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Appendix 1 

Integrable Potential Perturbations of the Ellipti¬ 
cal Billiard 


The equation 

>yVxy + 3 (j/t4 - XVy) + - x'^)V^y + xy (Kx - Vyy) = 0, (lO) 


is a special case of the Bertrand-Darboux equation [illliss], which represents 
the necessary and sufficient condition for a natural mechanical system with two 
degrees of freedom 

H = \{pl+pl) + V{x,y) 

to be separable in elliptical coordinates or some of their degenerations. 

Solutions of Equation (uni) in the form of Laurent polynomials in x,y were 
described in [TB]. Such solutions are in [TH] naturally related to the well-known 
hypergeometric functions of Appell. This relation automatically gives a wider 
class of solutions of Equation m - new potentials are obtained for non-integer 
parameters, giving a huge family of integrable billiards within an ellipse with 
potentials. Similar formulae for potential perturbations for the Jacobi problem 
for geodesics on an ellipsoid from min] are given. They show the existence 
of a connection between separability of classical systems on one hand, and the 
theory of hypergeometric functions on the other one, which is still not completely 
understood. Basic references for the Appell functions are [2l [3l I35j . 

The function F 4 is one of the four hypergeometric functions in two variables, 
which are introduced by Appell mm- 


F 4 {a,b, c, d-,x,y) 


(^)m-t-n(^)m-t-n X y 

^ {c)m{d)n ml nl ’ 


where (a)„ is the standard Pochhammer symbol: 


(a)n 


r(a -I- n) 
r(a) 


(a)o = 1. 


a(a + 1 )... (a + n — 1 ), 


(For example ml = (l)m-) 

The series is convergent for ^/x + ^/y < 1. The functions F 4 can be analyti¬ 
cally continued to the solutions of the equations: 


x(\ — x) 


y(i - y) 


d'^F 

dx"^ 


,d^F 


y 


dy-^ 


- 2 xy- 


d^F 


dxdy 

dF dF 

-I- [c — (a -I- 6 -I- l)a:] -- {a + b+ 1)?/-^- abF = 0, 

ox oy 


d^F ^d^F 


dy^ 


dx^ 


- ‘2‘Xy- 


d^F 


- + 


dxdy 

dF dF 

+ [c' — {a+ b+ l)y]-r-(a -I- & -I- l)a:—- abF = 0, 

dy dx 
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Al.l Potential Perturbations of a Billiard inside an Ellipse 

A billiard system which describes a particle moving freely within the ellipse 

2 2 

_ + ^ = 1 
A B 

is completely integrable and it has an additional integral 

f {xy-yxf 
+ AB • 

We are interested now in potential perturbations V = V(x,y) such that the 
perturbed system has an integral Ki of the form 

Ki = Ki + ki(x,y), 

where fci = ki{x,y) depends only on coordinates. This specific condition leads 
to Equation (fTTill on V with X = A — B. 

In [161117| the Laurent polynomial solutions of Equation CHD were given. 
Denoting 

V = -7)iEi(l,2 -7,2,1 - 7,i,y) + l), (11) 

where x = x^/A, y = —y^/A, the more general result was obtained in [18] : 

Theorem 11 Every function Ey given with (HD and 7 G C is a solution of 
Equation m- 

This theorem gives new potentials for non-integer 7 . For integer 7 , one 
obtains the Laurent solutions. 

Mechanical Interpretation. With 7 G ]R“ and the coefficient multiplying 
Ey positive, we have a potential barrier along x-axis. We can consider billiard 
motion in the upper half-plane. Then we can assume that a cut is done along 
negative part of y-axis, in order to get a unique-valued real function as a po¬ 
tential. 

Solutions of Equation m are also connected with interesting geometric 
subjects. 


A1.2 The Jacobi Problem for Geodesics on an Ellipsoid 

The Jacobi problem for the geodesics on an ellipsoid 


has an additional integral 


x^ 

~A 


B 


C 


1 


Ki = 


^ c^j 
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Potential perturbations V = V(x,y,z), such that perturbed systems have inte¬ 
grals of the form 


Ki = Ki + k(x,y,z), 


satisfy the following system: 


i + li + ilW. 




2 


(J2 

,2 


xy 


A- B ^ y X Vy 

AB B^ A A^ b' 


xy 


— -^]V + 

A3 B3 i + AB 


K 


vv 

A 


1 /, 


XX 

W 


CA^ 


V.y- 


zy 

CB^ 


K,. = 0 




A 2 


52 

2 




B-C Z Vy y V, 

BC B B^ C 




53 6*3 




yz 

BC 


14 


ZZ 

'W 


^yy 

C 


xy 

AB^ 


P..- 




( 12 ) 


x^ 2 /^ z'^ C — A „ X Vz 

_i_ ^ _I_I y _a_i 

A^ B^ Cy """ AC A2 C 


, z 14 

5- - 

C2 A 

zy 


(z^ x^\^^ xz f Vxx 14 z^ zy yx 

+ ^ ^ V~ ~ 

This system is an analogue of the Bertrand-Darboux equation m (see [) 
Let 


x‘^C{A-C) „ y'^C{C-B) 

z'^{B-A)A z^{B-A)B 

The following statement was also proved in [18]: 


Theorem 12 For every 7 S C, the function 

z'^V 

-F4(l,-7 + 2,2,-7-hl,i,y) 

is a solution of the system m- 

Thus, by solving Bertrand-Darboux equation and its generalizations, as it 
is done in Theorems ttH and m one get large families of separable mechanical 
systems with two degrees of freedom. It is well known that separable systems 
of two degrees of freedom are necessarily of the Liouville type, see [55] . 

Now the natural question of Poncelet type theorem describing periodic solu¬ 
tions of such perturbed billiard systems arises. It appears that again Darboux 
studied such a question, since in m, he analyzed generalizations of Poncelet 
theorem in the case of the Liouville surfaces. 


14 = (-7 + 1) 
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Appendix 2 

Poncelet Theorem on Liouville Surfaces 

In this section, we are going to give the presentation and comments to the 
Darboux results on the generalization of Poncelet theorem to Liouville surfaces. 

A2.1 Liouville Surfaces and Families of Geodesic Conics 

In this subsection, following we are going to define geodesic conics on 
an arbitrary surface, derive some important properties of theirs and finally to 
obtain an important characterization of Liouville surfaces via families of geodesic 
conics. 

Let Cl and C 2 be two fixed curves on a given surface S. Geodesic ellipses 
and hyperbolae on S are curves given by the equations: 

9 + a = const, 

9 — a = const, 

where 9, a are geodesic distances from Ci, C 2 respectively. 

A coordinate system composed of geodesic ellipses and hyperbolae joined to 
two fixed curves is orthogonal. In the following proposition, we are going to 
describe all orthogonal coordinate systems with coordinate curves that can be 
regarded as a family of geodesic ellipses and hyperbolae. 

Proposition 8 Let 

ds"^ = A^dv^ + C^dv"^ (13) 

be the surface element corresponding to an orthogonal system of coordinate 
curves. Then the coordinate curves represent a family of geodesic ellipses and 
hyperbolae if and only if the coefficients A, C satisfy a relation of the form: 

^2 + c-^ - ’ 

with U and V being functions of u, v respectively. 

Proof. By assumption, equations of coordinate curves are: 

9 + a = const, 9 — a = const, 

with 0 , a representing geodesic distances from a point of the surface to two fixed 
curves. 

Thus: 

u = F{9 + (t), V = Fi{9 — a). 

Solving these equations with respect to 9 and a, we obtain: 

9 = 4>{u) + il)(v)., a = <f)(u) — ifiv). 


43 


As geodesic distances, 9 and a need to satisfy the characteristic partial 
differential equation: 

^ \dv> + V 3 m 1 ('14', 

A2C2 ^ ’ 

From there, we deduce the desired relation with U = {(j)'{u))^, V = {ip'{v))‘^. 
The converse is proved in a similar manner. □ 

As a straightforward consequence, the following interesting property is ob¬ 
tained: 

Corollary 3 If an orthogonal system of curves can be regarded as a system of 
geodesic ellipses and hyperbolae in two different ways, then it can be regarded as 
such a system in infinitely many ways. 

Now, let us concentrate to Liouville surfaces, i.e. to surfaces with the surface 
element of the form 


ds^ = {U-V){Uidu^+ V^dv^), (15) 

where U, Ui and V, Vi depend only on u and v respectively. 

Now, we are ready to present the characterization of Liouville surfaces via 
geodesic conics. 


Theorem 13 An orthogonal system on a surface can be regarded in two differ¬ 
ent manners as a system of geodesic conics if and only if it is of the Liouville 
form. 


Proof. Consider a surface with the element (fT^ . If coordinate lines u = const, 
V = const can be regarded as geodesic conics in two different manners then, by 
Proposition 6, A, C satisfy two different equations: 


— — - 1 
^2 + C'2 - ’ 


El 

A2 



Solving them with respect to A, C, we obtain that the surface element is of the 
form 

Conversely, consider the Liouville surface with the element dni) and the 
following solutions of the equation (fHl) : 


9 = 

a = 


Ui\/U — a du-\- j V\\Ja — Vdv, 
UiVU — adu — J ViVa — Vdv. 


The equations 9-\-a = const, 9 — a = const will define the coordinate curves. □ 
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A2.2 Generalization of Graves and Poncelet Theorems to 
Liouville Snrfaces 

We learned from Darboux m that Liouville surfaces are exactly these having 
an orthogonal system of curves that can be regarded in two or, equivalently, 
infinitely many different ways, as geodesic conics. Now, we are going to present 
how to make a choice, among these infinitely many presentations, of the most 
convenient one, which will enable us to show the generalizations of theorems of 
Graves and Poncelet. All these ideas of enlightening beauty and profoundness 
belong to Darboux m 

Consider a curve 7 on a surface S. The involute of 7 with respect to a point 
A G 7 is the set of endpoints M of all geodesic segments TM , such that: 

T 

TM IS tangent to 7 at T; 

the length of TM is equal to the length of the segment TA C and 

these two segments are placed at the same side of the point T. 



Figure 13: Involute 


Involutes have the following important property, which follows immediately 
from the definition: 

Lemma 10 The geodesic segments TM are orthogonal to the involute, and the 
involute itself is orthogonal to 7 at A. 

Now, we are going to find explicitly the equations of involutes of coordinate 
curves on a Liouville surface S with the surface element (Ha. 

Lemma 11 The curves on S given by the equations: 


^ — J UiVU — a du + j Vi^/a — Vdv = const, 
a = J Ui \/U — a du — J Vi Va — Vdv = const 


(16) 


are involutes of the coordinate curve whose parameter satisfies the equation: 

[U -a){V -a) =0. 
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Proof. Fix the parameter a. The equations of of geodesics normal to the curves 
9 = const, a = const are obtained by differentiating m with respect to a: 


Uidu Vidv 
y/U — a y/a — V 


(17) 


Let uq be a solution of the equation U — a = 0. Then, the geodesic line inn 
will satisfy du = 0 at the point of intersection with the curve u = uq, i.e. it will 
be tangent to this coordinate curve. The statement now follows from Lemma 

[ini □ 


Proposition 9 Coordinate curves on a Liouville surface are geodesic conics 
with respect to any two involutes of one of them. 

Proof. Follows from Lemma [TT] and the proof of Theorem [131 D 

Now, we are ready to prove the generalization of Graves’ theorem. 

Theorem 14 Let £q : u = uq and £i : u = ui be coordinate curves on the 
Liouville surface S. For a point M G £i, denote by MP and MP' geodesic 
segments that touch £q at Q, Q'. Then the expression 

e{MP) + £(MP') - £{PP') 

is constant for all M, where i{MP), £{MP'), and £{PP') denote lengths of 
geodesic segments MP, MP', and of the segment PP' C £q respectively. 



Figure 14: Theorem HM 


Proof. Let T>, V be involutes of the curve £o with respect to points R, R' G £q, 
and Q,Q' intersections of geodesics MP,MP' with these involutes. 

Both £o and £i are geodesic ellipses with base curves T), V, thus the sum 
£{MQ) + £{MQ') remains constant when M moves on £i. 

Since £{PR) = £{MP) + £{MQ), £{P'R') = £{MP') + £{MQ’), we have: 

£{MQ) + £{MQ’) = £{PR) + £{P'R’) - £{MP) - £{MP') 

= £{RR') - {£{MP) + £{MP') - £{PP ')), 


46 









and the theorem is proved. 


□ 


From here, the complete analogue of the Poncelet theorem can be derived: 

Theorem 15 Let us consider a polygon on the Liouville surface S, with all 
sides being geodesics tangent to a given coordinate curve, and each vertex but 
one moving on a coordinate curve. Then the last vertex also remains on a fixed 
coordinate curve. 
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